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Abstract. We present a new method to parameterize Type la Supernovae (SN la) multi-color light curves. The method was 
developed in order to analyze the large number of SN la multi-color light curves measured in current high-redshift projects. 
The technique is based on empirically modeling SN la luminosity variations as a function of phase, wavelength, a shape 
parameter, and a color parameter. The model is trained with a sample of well-measured nearby SN la and then tested with an 
independent set of supernovae by building an optimal luminosity distance estimator that combines the supernova rest-frame 
luminosity, shape parameter, and color reconstructed with the model. The distances we measure using B- and V-band data show 
a dispersion around the Hubble line comparable or lower than obtained with other methods. With this model, we are able to 
measure distances using U- and S-band data with a dispersion of 0.16 ± 0.05 around the Hubble line. 

Key words, supernovae: general - cosmology:observations 



1. Introduction 

Type la supernovae (SNe la) are a powerful tool for studying 
the evolution of the luminosity distance as a function of red- 
shift and for subsequently constraining the cosmological pa- 
rameters. SNe la are indeed very luminous and "standardiz- 
able" candles, a nd have l ead to the discovery of the accele ration 
of the Universe jRiess et alll998tlPerlmutter et alJl999l) . 

Although often described as a homogeneous class of ob- 
jects, SNe la exhibit variability in light curve shapes, colors, 
intrinsic luminosity, and spectral features. Finding correlations 
among SN la observables is motivated by improving the esti- 
mation of their intrinsic luminosity on an event-by-event basis, 
in order to reduce the scatter in luminosity distance estimates. 
The main correlations observed in photometric measurements 
are: 

- a width-luminosity (or brighter-slower) relation, which ex- 
presses the fact that brighter supe rnovae have a slower de- 
cline rate than fainter on es (Pskovskii 1977; Phillies 1993; 
Riess et alJ Il995t lHamuv et alJ Il996bl IPerlmutter et alJ 
19971) . 

- a brigh t er-bluer relation which was made exp l icit inlTripd 
Jl998l) : iTripp & Branch! J 19991) : IParodi et all feOOOh . and 
assume d to be due to extinct i on by dust in other 
works JRiess et alJ Il996al 1 19981: IPerlmutter et alJ 1 1999; 
iTonrv et all200.*lKnor> et alfcOOSHBarris et al l2004). 

The case of the brighter-bluer relation is interesting. Even if 
authors fully agree neither on the origin of the effect nor on the 
strength of the correlation, most, if not all, recent attempts to 



build a SN la Hubb le diagram have ma de use of color in their 
distance estimator 1 . iRiess et alJ il996b) summarizes previous 
work on the subject and proposes a way to reconcile diver- 
gent interpretations of data by taking into account the corre- 
lation between light curve shape and color. New methods us- 
ing color information have also been recentl y proposed to es- 
tima te luminosity dista nces (see for example Wan g et al.l f2003 
and lWang etal]|2005i) . We will bring to the debate our own 
estimations of the correlation strength and its wavelength de- 
pendence. 

Cosmological measurements using SNe la are based on 
comparing nearby and distant objects. In order to reduce the 
sources of systematic uncertainties, it is important that all dis- 
tances are derived using the same procedure, especially when 
considering large samples such as those being collected by the 
ESSENCE 2 or SNLS 3 projects. To analyze these data sets, the 
following constraints have to be taken into account: 

- High redshift objects often lack late-time photometric data 
(or have one of too poor quality) which makes it impracti- 
cable Jx^stimate color from late-time data as proposed by 
lLiraetalJ ( ll998l) . 

- Above redshifts z ~ 0.8, rest-frame [/-band measurements 
have to be used because of the limitation of silicon detec- 
tors. Incorporating rest-frame U -band measurements on the 



1 In lPerlmutter et all 1 1999). it is checked that B — V rest-frame col- 
ors of nearby and distant type agree on average, and the color mea- 
surement is not used event per event. 

2 http://www.ctio.noao.edu/~wsne 

3 |http : //c fht . hawaii . edu/S NLS| 
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same footing as the commonly used rest-frame B- and V- 
band data is hence highly desirable. 
- Applying cuts on light curve parameters should be prohib- 
ited since the resolution on these parameters follows the 
photometric resolutions and hence degrades with redshift, 
therefore biasing the event sample. In particular, cuts on 
colors to eliminate reddened events are to be avoided. 

Various tec hniq ues have been used to e stimates distances. 
IPhillipsI dl993l) and IPerlmutter et all dl997l) fit one band at a 
time, and derive distances from light curve parameters mea- 
sured in different bands (usually B and V bands), sometimes 
relying on late time measurements to measure host galaxy ex- 
tinction. Are^ejiUmdrejjried vers ion of this approac h can be 
found in Wa ng etalJ J2005I) . and in lWang et aD 00031) it is pro- 
posed to fit B - V as a function of B. These m ethods do not 
make use of rest-frame £/-band measurements. iNugent et all 
discuss in detail the problems associated with SN la 
[/-band photometry. They mention in particular the apparent 
large intrinsic variations of the UV luminosity among similar 
supernovae as well as our lack of understanding of the SN la 
UV photometry -principally due the uncertainties which affect 
the large extinction corrections which have to be applied to the 
data. Nevertheless, it is important to try and incorporate rest- 
frame [/-band data t o estimate distances. 

An extension of Perlm utter et alJ dl997h to rest-frame U- 
band data was developed in Kno p et al.l (2003), but at the ex- 
pense of adding a large distance systematic (and probabl y sta- 
tistical) unce rtainty. Similarly, the MLCS method dRiess et alJ 
Il995t l 1996a) has been extended to include [/-band measure- 
ments under the n ame of MLCS2k2 (Jha 2002) and used in 
iRiess et all d2004l). but also at the expense of a worsened dis- 
tance resolution ( Jha 120021) . We will show that our model pre- 
dicts light curves for any band located between rest-frame U 
and R bands and that we are able to get a similar or better dis- 
tance resolution with both rest-frame (B, V) and ( U,B) band 
pairs. 

We propose here to parameterize the light curve model with 
a minimal parameter set: a luminosity parameter, a decline rate 
parameter and a single color parameter. Our approach will be 
to build a phenomenological model of the expected SN flux, 
continuously varying with phase, wavelength, decline rate and 
color, in order to capture all these features at once. This ap- 
proach offers several practical advantages which make it easily 
applicable to high-redshift SNe la currently measured in large 
projects. First, the k-corrections are built into the model and 
not applied to the data. This allows one to propagate all the un- 
certainties directly from the measurement errors. More impor- 
tantly, when needed, we make use of the SN rest-frame U -band 
fluxes to estimate the supernova distances. 

In Sect. 2 we describe the semi-analytic model used. We 
then describe, in Sect. 3, how the coefficients of the model 
are determined by an iterative training based on a set of well- 
sampled nearby SNe la taken from the literature. We also high- 
light some properties of the resulting model. At this stage, the 
aim is to model multi-color light curves and not to estimate 
luminosity distances. The model is tested in Sect. 4 with an in- 
dependent set of SNe la in the Hubble flow. A luminosity dis- 



tance estimate is then constructed from the fitted parameters of 
the light curve model. It is used to build Hubble diagrams suc- 
cessively from (B, V) and (U, B) light curve pairs. In order to 
assess the precision of this approach, we compare distance es- 
timates of the same events obtained from the (U, B) and (B, V) 
band pairs and with other distance estimators. 

2. The light curve model 

2.1. Model definitions 

As already mentioned, we choose to parameterize light curves 
(more precisely light curve pairs or light curve triplets when 
available) using a single luminosity, a single shape parameter 
and a single color. The choice among possible implementations 
is largely arbitrary. We choose the following parameters, which 
enable comparisons to be made with previous works: 

- f : a global intensity parameter which varies with redshift 
like the inverse of the luminosity distance squared, 

- s : a time stretch fa cto r as the decline rate indicator 
dPerlmutter et al.ll997t) . In^pldhaber etaPd200ll) . this pa- 
rameter is shown to apply to the rising part of the light curve 
as well. However, while the stretch paradigm describes well 
the bright part of the B light curve, it does poorly at late 
time. It also fails to capture the shape variations in the other 
bands. This is why our model uses the stretch parameter 
as an index rather than the stretch paradigm itself. As de- 
scribed below, by construction our model follows exactly 
the stretch paradigm in the B band. 

- c = (B - V) max + 0.057, where (B - V) max is measured at B 
maximum, and -0.057 is the chosen reference color (Vega 
magnitudes) of a SN la. 

- t^ wx : the date of maximum in the rest-frame B band 

With these definitions, the expected counting rate /sn in a given 
pass-band T, of a supernova at redshift z, and at a phase p = 
(t - tf nax )l(\ + z), can be written: 

f SN (p, Z, T)=f (l+ z) J" 4>(p, A, s, c) £ T(A(1 + 2)) dA 

<p{p, A, s, c) is a model of the SN la energy luminosity per unit 
wavelength. It may vary with the supernova stretch and color. 
Note that the potential extinction by dust in the host galaxy is 
not explicit in the equation. Instead, we choose to incorporate 
it in the model (p{p, A, s, c) as discussed below. 

Building an average spectral template as a function of 
phase, wavelength, color and stretch from observations is com- 
plicated because of the inhomogeneity and incompleteness of 
published data. Although some spectral features have been cor- 
related with stretch, we do not have yet a complete knowledge 
of spectral diversity as a function of phase. In most of the ap- 
proaches to SN la light curve fitting, the limited knowledge 
of spectral variability impacts on the accuracy of cross-filter 
k-corrections, defined as expected ratios of fluxes in differen t 
bands at the same phase ijKim et alJl996HNugent et all2002l) . 

In lNugent et all d2002T) T it was shown that the variation in 
k-corrections from one supernova to another depends primar- 
ily on the supernova color, and to a lesser extent on spectral 
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features 4 . Hence we neglected the variability of those spectral 
features in the modeling of light curves. 

In order to implement stretch dependent light curve shapes 
and colors we therefore used the following approximation: 

fs N (Ps,Z, r) = / (l+z) J" <p(p s ,A)-£-T(A(l+z))dA 

x exp [-0.4 ln(10) x 9C(p s , A T , s, c)] (1) 

where p s = p/s is a stretch-corrected phase. This functional 
form defines the light curve model. 

In equation^ <t> no longer depends explicitly on s and c, and 
< K{p s , A, s, c) is a smooth "correction" function of our four vari- 
ables. At is the central wavelength of the filter T. 'K enables one 
to implement light curve shape variations that are more compli- 
cated than simple dilation of the time scale, along with stretch 
dependent colors. As described below, 'K varies smoothly with 
A; this justifies placing it outside the integral over wavelength. 
We also considered keeping "Kips, At, s, c) inside the integral. 
It changed the results of the fits by less than 1%, while the 
computing time was multiplied by a factor of 10. 

Equation Q implements k-corrections for an average 
SN la, conforming to the common practice. The approach pro- 
posed here will be easy to adapt when constructing a stretch- 
dependent spectral template becomes possible. 

The k-corrections are usually applied to data. Here, they 
are incorporated into the model. This offers a few practical ad- 
vantages: first, light curves can be generated for arbitrary pass- 
bands (within the spectral coverage of <\> and "K); light curves 
in the observed pass-bands can be directly fitted to the data. 
Second, the light curve parameter uncertainties extracted from 
the fit incorporate all uncertainties propagated from measure- 
ment uncertainties; for example, uncertainties introduced in the 
k-corrections by the possibly poor determination of the date of 
maximum and/or color are propagated into the parameters un- 
certainties. 

The functions (f> and 'K define the model. Once they are 
determined, one can fit the supernova photometric data points, 
measured in a minimum of two pass-bands, to estimate fo, s, c, 
and a date of B maximum light, which is a nuisance parameter. 
With only one passband, c must be held fixed. 

For (pjp ., A), we use a tem plate spectrum assembled by 
P. Nugent ( iNugent et alJ l2002l and private communication) 
smoothed along the phase (time) axis, and normalized as a 
function of phase to the g-ban d light curve template "Parab 
-18" of Goldhaber et alJ 11200 II) . Any smooth variation of the 
template (f> with phase or wavelength is irrelevant at this level, 
since it may be changed by the function 'Kips, A, s, c). The only 
important quantities in (p(p s , A) are the SN la spectral features, 
intended to be realistic on average. 

The empirical correction function 7C is implemented as a 
sum of two polynomials: 

K(p s , A, s, c) = K s {p s , s) + n< c {A, c) (2) 

4 Most of the SN la photometric reductions transform instrumental 
magnitudes into standard magnitudes using color equations, derived 
from standard stars observations. This assumes that color rather than 
spectral features dominates cross-filter corrections. 



where we explicitly separate the corrections associated with the 
parameters s and c to clarify their interpretation. ( K s (p s , A, s) 
modifies the shape of light curves and absorbs any stretch- 
color relation except for the (B - V) max color. Indeed we want c 
to describe exactly the (B - V) max color. *7C C (A, c) is then a color 
correction as a function of wavelength and color. 

In order to remove all degeneracies among coefficients, the 
model must fulfill the following constraints: 

9C s (p s < 35, s,A B ) = 

9C s (p s = 0, s, A v ) = 

<K c (A,c = 0)=<K c (A B ,c) = 

<K c (A B ,c)-<K c (A v ,c) = c 

where A B and Av refer to the mean wavelengths of the B and V 
filters. The first constraint ensures that the parameter s actually 
defines the stretch in B band, since the template is not modified 
for A = A B . The remaining constraints define c as the B - V 
color (relative to the spectral template) at maximum B, for all 
stretches. Note that other colors may (and in fact will) depend 
on s at fixed c. fit describes the actual peak flux in B. It is in no 
way "corrected" for the brighter-slower relation. 

7f s is implemented as a polynomial of degree D p in phase, 
Dx in A and D s in stretch respectively. Similarly, % c is a poly- 
nomial of degree D p in phase and D c in color. In this study, 
we chose for the degrees of polynomials (D p , Dx, D s , D c ) = 
(4, 3, 1, 1). The large degree in phase permits a detailed adap- 
tation of the model light curve shapes to the actual data. The 
other degrees correspond to the minimal number of coefficients 
(a degree 3 in wavelength is chosen to adapt colors in UBVR, 
independently). The % polynomials are then defined by 48 co- 
efficients, but the constraints reduce the number of independent 
coefficients to 34. There are 32 free coefficients for % s and 2 
free coefficients for 7C C . We will call "training" the determina- 
tion of these coefficients from measurements of nearby SNe la. 

2.2. Normalization of transmissions 

The value of fo depends on the normalization of the spectral 
template (f>. This is not an issue when comparing nearby and 
distant supernovae to measure cosmological parameters, since 
only flux ratios matter: all objects have to be fitted using the 
same model. For Hq measurements, one would need to nor- 
malize (0) using SNe la at known distances, but we will not 
perform this here. 

More important is the apparent dependence of the value 
of fo with respect to the normalization of the transmission T 
(equation^. Since we aim at fitting multi-color light curves, 
observed with different instruments, with a single set of param- 
eters (fo, s, c), it is mandatory to remove this dependence. 

All photometric data are expressed in units of the integrated 
flux f re j (deduced from the zero point) of a known standard 
spectrum <p, e f (e.g. Vega). While the functional form of T(A) is 
determined by optical transmission measurements its normal- 
ization can be determined from f re f and <p re f via the relation 

J <p ref (A)-?-T(A)dA = f ref (3) 
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The filter transmissions are usually published as series of 
numbers (as inlBessell Il990l). and this leads to ambigu ities: 
for ex ample, ISuntzeff et al.l d 19991) writes "Note that IB esse lA 
\l99(k defines the sensitivity function as the product of the 
quantum efficiency of the detector+telescope, the filter trans- 
mission curve, the atmospheric extinction, and a linearly in- 
creasing function of wavelength." . The ambiguity is whether 
the transmission refers to signal per unit energy, or signal per 
photon flux. Attributing dimensions to transmissions solves the 
ambiguity. 

There have also been concerns about the relative weights 
of the different wav elengths in the integrals over wavelength 
(Nug ent et al.l2002l) . namely whether one should sum photons 
or energies. If one aims at reproducing the instrument response 
to an arbitrary spectral energy density, the mathematical inte- 
grals sho uld mimic the p hysical integration process of the in- 
strument (Fukugitoetaljl993). For example, if one considers a 
CCD-based observing system, the effective transmission T(A) 
will be proportional to a number of photo-electrons per pho- 
ton. So one should not integrate photon counts nor energies, 
but charge on the detector (or ADC counts). 

3. Training the model 

Since /o describes the observed luminosity in the B band (be- 
cause of the constraints applied to the corrections), the model 
only incorporates stretch-shape and stretch-color relations, but 
no correlation involving luminosity. This option was chosen in 
order to allow us to train the model with objects at unkn own 
distanc es, in particular the nearby objects in the sample of l.Thal 
(2002) measured in the U band. If one offsets all magnitudes 
of each training object by an arbitrary amount, possibly differ- 
ent for each object, the resulting model will not change. One 
could then consider incorporating high redshift objects into the 
training, but we choose not to do it here. 

3.1. The nearby SN la Sample 

The model was trained and tested using a sample of published 

nearby supernova light curves. We collected 122 SNe la for 

which B- and V-ban d light curves are available i n the li t erature , 

I n 1 1 — - — I I 1 j — - — i 

including data from Hamuv et al. ( 1996b); Riess et al. ( 1996b) 

and Jha (2002) for a total of 94 objects, and 28 additional su- 
pernovae collected from various sources (see the caption of ta- 
ble[l). 

Objects were then selected based on two main criteria. 
First, we kept supernovae with at least two measurements be- 
fore the maximum in the B or the V band. This is necessary 
to ensure that the date of maximum is well defined and that 
the measurements can safely be used as a function of phase. 
Out of the whole sample, 56 SNe satisfied this criterion. Then 
under-luminous peculiar sup ernovae: SN 1991bg, SN 1998bp, 
SN 1998de and SN 1999bv (lHowelll200lUT,i et al feOOIhl). and 
the peculiar objects SN 2000cx and SN 2002eri lLi et al 12001 at 
Pign ata et alj f2004) were rejected from our sample. There are 
a number of reasons for this choice; the most important one 
being the spectral difference between normal and SN 1991bg- 
like SNe la. Since our model is built using a spectral tem- 



plate describing the features of the average normal SN la, it 
is not well suited for describing very different objects, such as 
SN 1991bg-like events. We note that, this is not a constraint 
for cosmological studies since, so far, no under-luminous SNe 
have been found in distant SN searches. Furthermore, if present 
at high-redshift, these objects would easily be identified both 
from spectr oscopy and from their B - V color evolution (see 
for instance Iphillips et all 1 9991 iGarnavich et all2004T) . 

SN 1991T-like events were kept in the sample. This "sub- 
class" of over-luminous SNe la is spectroscopically identified 
by the presence of unusually weak absorption lines during 
the pre-maximum phase. These are more difficult to identify 
than their under-luminous counterparts in low signal-to-noise 
spectra, such as those usually available for high redshift SNe. 
Moreover, the B - V color evolution is not particularly different 
from the one of normal SNe. Ideally, one should use a different 
spectral template for fitting this kind of events. We estimate the 
errors introduced in the k-corrections using a standard SN la 
spectral template to be of a few percent in the worst cases, 
and well within the final dispersion in the peak luminosity of 
SNe la. We further note that, the uncertainty due to using a non 
optimal spectral template is analogous to the uncertainty in the 
cross filter ^-corrections used in standard methods. 

The resulting sample of 50 SNe was then split into two sets: 
a training sample and a test sample. The training sample was 
used to adjust the coefficients of the polynomials of the model. 
It contains all the supernovae with redshifts smaller than 0.015 
(not in the Hubble flow) and 6 supernovae at redshifts above 
0.015, for which with [/-band data was available, in order to 
improve the model in this wavelength region. The training sam- 
ple contains the 34 supernovae listed in tabled The test sample 
contains 26 supernovae (table|2j. Note that the two samples are 
not completely independent. Indeed, they share ten supernovae 
with U -band light curves, since such events are rather scarce. 

The data was not pre processed in any way prior to fit- 
ting. To account for the Milky Way extinction, we incorpo- 
rate it into the ins trument transmission, using the law from 
ICardelli et al. (1989) with a color excess E(B — V) obtained 
from lSchlegel et al.ldl998t) dust maps at the position of the ob- 
ject to fit. 

3.2. Training the Model 

All the published nearby supernova magnitudes are expressed 
in the Johnson-Cousins UBVR system. In Equation Q, we 
use models of the instrument transmissions as a function of 
wavelength. We a dopted the transmission functions published 
bv lRessellNl990ll . and interpreted them as A T(A) (see equa- 
tion Q}, i.e. counts p er unit energy, following a footnote of 
ISuntzeff etaH l ll999l) . 

Training the model consists in determining the 
Kip.?, A, s, c) = < K s (p s , A, s) + ftci^, c) correction func- 
tion (Eq.|3 using the training sample data in the UBVR bands. 
The R band was introduced to avoid relying on extrapolation 
for measurements beyond rest-frame V band. We start with a 
first guess: 

%(p s ,A,s) = 
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K c (A,c) = cx(A-A B )/(A v -A B ) 

and we use an iterative algorithm which can be sketched as 
follows: 

1 . Fit the light curves using the current determination of 'K. 

2. Fit 6 ( K(p s ,A,s,c), an instance of the "7C function, on the 
light curve residuals. During this step, identify and remove 
the outliers data points. 

3. K + 5K. 

4. GOTO step 1, until 6 7C becomes negligible. 

All the available data, i.e. the UBVR residuals, up to 
large phases were used to determine the 7C function (step 2). 
However, not all the data points were used to fit the light curves 
(step 1): the /?-band data was not used and only U BV photomet- 
ric points with phases ranging from -15 days to +35 days were 
used since we are mostly interested in describing the rest-frame 
UBV central part of the supernova light curves. 

3.3. Results of the training 

The fit converged after four iterations. 2480 measurement 
points were fitted, and 39 were discarded as outliers (at the 
3 cr level). Compared to the number of free coefficients of 
the model, we can safely conclude that the model is not over- 
trained. The standard deviations of the residuals to the model 
in UBVR are respectively of 0.09, 0.09, 0.06, 0.07 magnitudes. 
Figure Q shows the final U, B, V and R templates obtained at 
the end of the process as a function of stretch. By construction, 
the rest-frame B and V-band magnitudes at maximum do not 
vary with stretch. We find a strong dependence of (£7 — B) m „ x 
with stretch (6(U - B) max -6s, compatible with ljhal Eb02). 
which is an essential feature for the model to reproduce in or- 
der to estimate a reliable color in the wavelength range between 
U and B. The model also manages to reproduce a a stretch- 
dependent secondary shoulder in the R band. We also notice 
that the /?-band light curves cross each-other for different val- 
ues of the stretch. This reproduces well the fact that the brighter 
slower relation i s weak er in the redder pass-bands, as noted 
bv lPhillips et alJ dl999). The residuals to the light curve fit are 
shown figure |2 

Figure[3]represents the color correctio n < K r (A, c) for c = 0. 1 
compared to the dust extinction law from Card elli et alJ (fl989). 
Interestingly enough, the law we obtain follows pretty well that 
of Cardelli in the R band but not in the U band where we get 
a stronger dependence on c. Also shown (shaded area) is the 
uncertainty on 7C C derived from the x 1 increment (normalized 
to the number of supernovae) of the fit to the light curve residu- 
als. The Cardelli law in the U band is at 3.7 standard deviation 
from the best fit value. 

As a consequence, we deduce that the relation between 
E(B-V) and E(V-R) are very similar (i.e. indistinguishable) 
to the o nes expected from r eddening by dust. This similarity, 
noted bv lRiess et alJ(ll99 6b). does not prove however that c can 
be interpreted as reddening by dust; an additional requirement 
for this hypothesis to be valid would be that the peak Z?-band 
magnitude increases with c by a value of Rb x c. We will see 
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Fig. 1 The UBVR template light curves obtained after the train- 
ing phase for different values of the stretch and null color. Note 
the strong variations of the f/-band maximum with the stretch. 
Note also the after max shoulders in the V and R bands. The 
model also reproduces well the fact that the brighter-slower re- 
lation is weaker in the R band than in the B or V bands. The 
model does not incorporate the brighter-slower correlation. The 
apparent brighter-slower correlation for the U band is, in fact, 
a stretch-color correlation. 



that this is not the case in the next section. Let us also empha- 
size that the stretch dependent part of the U-B and V-R colors 
are included in the stretch dependent term 'KsiPs, A, s), and not 
in the color curve of figure[3] 



4. Performance study 

Once the correction function 'Kips, A, s, c) is determined, we 
can fit the model on the sample of nearby SN la light curves 
listed in table |2] This allows us to perform various consistency 
checks, in order to make sure that the model describes well the 
UBV photometry of SNe la. 
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10 20 
phase 



Fig. 2 Residuals to the light curve fit as a function of phase, 
for supernovae from the training sample (open symbols) and 
the test sample (filled symbols). There is no significant dif- 
ference between the distributions obtained with both samples 
which proves that the model is not over-trained. The systematic 
residuals as a function of phase are a direct consequence of the 
limited number of parameters (5) used to implement the correc- 
tions TO, as a function of phase. Increasing this number would 
remove this systematic trend in residuals, but some parameters 
would be poorly constrained in the U band where photometric 
data is scarce. 



A first test consists in checking the ability of the model to 
reproduce the shape and color features of the independent set 
of SNe la. This is demonstrated on figure |2 which shows the 
residuals to the fits of light curves of SNe from the test sample. 



4.1. Distance estimate 

The global intensity parameter /o is proportional to djXzf 2 and 
inversely proportional to the normalization of 0. One can define 




4000 



5000 

MA) 



6000 



Fig. 3 The color correction 7C e as a function of wavelength for 
a value of c of 0.1 (thick line). The shaded area corresponds to 
one standard deviation to the best fit value. The dashed curve 
repre sents the extinction with respect to B band, (A^ - A B ), 
from lCardelli et alJ i 19891) with R v = 3.1 and E(B-V) = 0.1. 
The curves are compatible in the R band but not in the U band 
(by construction they match in B and V bands as illustrated 
with the open circles). 



a rest-frame B magnitude m* B JPerlmutter et al.l 19971) which re 
moves this artificial dependence on the model normalization, 

fsN(0,Z,T* B ) 



m* = -2.5 log I0 ■ 



(4) 



(l+Z)fref(T B ) 

where and f re f are respectively defined by equations ^ 
and|3j Tg is the transmission of the B filter and 
T B (A) = Tb(A./(1 + z)) is a redshifted B transmission. One 
can check that m* B varies as 5 log 1() ddz) with redshift and that 
m* B —> nig for z <K 1, where ma is the conventional B magni- 
tude. 

We incorporate the Hubble parameter dependence of di in 
a constant parameter M™ = Mb - 5 log 1() (/i7o), which is the 
average absolute magnitude of a SN la with s — 1 and c — 0, 
for a value of the Hubble parameter of 70 km.s" 1 .Mpc~ I . 

As mentioned in the introduction, the peak luminosity of 
SNe la is correlated to stretch and color, so we may build a 
distance estimator that accounts for tho se correlation s and as a 
result reduces the dispersion. Following Tripp ( 1998), we adopt 
linear corrections of coefficients a and B respectively for stretch 
and color. The distance estimator is then 



-70 



i B ~M B 



- 43.16 + a (s-l)-Bc 



(5) 



Its expectation value for a supernova at redshift z is 
5 log 10 (cIl(z) Hqc~ 1 ^. Since our goal is here to test the distance 
estimator rather than actually perform a cosmological fit, we 
impose the "concordance" cosmological parameters (Om = 0.3 
and £2a = 0.7) when fitting Ml , a and B and apply the method 
to build low-z Hubble diagrams using successively (B, V) only 
and (U, B) only light curves of supernovae from the test sam- 
ple (table|2j- In general, rest-frame (U, B, V) light curve triplets 
should be used, when available, to determine the light curve 
parameters. 
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4.1 .1 . Hubble diagram in BV only 



Using B- and V-band only light curves of supernovae with 
redshifts larger than 0.015 from the test sample, we obtain: 
M™ = -19.41 ±0.04, a= 1.56 ±0.25 and /? = 2.19 + 0.33. The 
standard deviation of residuals is 0.16 ± 0.03 . Uncertainties 
on Mb, s, c along with their covariance were included in the fit 6 , 
we also considered an uncertainty on redshifts due to peculiar 
velocities of 300 km.s -1 ; an additional "intrinsic" dispersion of 
0.13 is needed in order to get a^- 2 per degree of freedom of 1. 
The resulting Hubble diagram as well as the observed brighter- 
slower and brighter-bluer relations are shown figure 0] 

Our approach t o estimating distances easily compares to 
the one adopted in lTrippI lll998h : the main differences are the 
light curve model and the brighter-slower parameterization. 
When w e fit the same SNe sample (The Calan-Tololo sam- 
ple from Ha muv et alJll996ah . a value of a — 1 .04 ± 0.24 and 
B = 2.08 ± 0.27 are obtained, which compare well to a = 0.88 7 , 
B = 2.09 of lTrippH 19981) . based on peak lumino sity, color, and 
decline rate estimates from lHamuv et al.(lfl99 6a). We conclude 
that our model correctly reproduces basic parameter estima- 
tions of previous works. 

Concerning the interpretation of the brighter-bluer correla- 
tion, we find a value of B which is incompatible with Rb = 4.1, 
expected for extinction by dust analogous to the observed 
law in the Milky Way. The value we find is c ompatible with 
those found in previous works (see Tripp ( 1998) and references 



therein). However, as stressed in lRiess et al. (1996b), the color 



excess (or deficit) at maximum should not be interpreted as 
entirely due to extinction but be corrected for the part of this 
excess that is correlated with stretch . We measure a stretch- 
color slope o f abou t 0.2, similar to the relation proposed in 
IPhillips et all dl999h 8 and can redefine our parameters to ac- 
count for this correlation: 

c = c + 0.2(5-1) 
s' = s 

so that s' and c' are uncorrected. The correlation coefficients 
then become a' — a + 0.2 B and R' = /3, which means that re- 
defining the color excess to explicitly assign to stretch the color 
variations correlated to stretch does not change the brighter- 
bluer correlation strength. 

The two-parameter (stretch and color) pragmatic approach 
we followed can accommodate both reddening by dust and any 
intrinsic color effect dependent or not on stretch. One may rea- 
sonably assume that reddening by dust and stretch ind ependent 
intrinsic colors mix (as proposed in lNobili et alJ2003i) . and that 

5 Note that this number takes into account the number of parameters 
in the fit. The measured RMS value is 0.14 ± 0.03. 

6 The uncertainties on the distance estimate formally depend on a 
and p, and increase with them. As a consequence, the x 2 minimum 
is biased toward large values of these parameters. We therefore com- 
puted the uncertainties with the initial values, and use the result of the 
fit at the final iteration. 

1 b = 0.52 translates to a ^ 0.88 when using stretch and the first 
order relation (AM 15 - 1.1) - 1.7(1 - s). 

8 The proposed relation is j^j- = 0.1 14 + 0.037. With the approx- 
imate relation 




0.05 
z 



0.1 





-1.7 (at s B = 1), we expect dc/ds ^ -0.2. 



Fig. 4 From top to bottom: Residuals to the Hubble diagram 
as a function of redshift, stretch s and color c indexes for su- 
pernovae of the test sample fitted in B and V bands. SNe with 
redshifts smaller than 0.015 are labeled with opened symbols. 



disentangling the contributions would improve the distance res- 
olution. Our distance indicator is however independent of the 
interpretation of the color variations. Since the low value of B 
may indicate that some intrinsic effect plays a role, we did not 
interpret color as only due to reddening by dust and hence ac- 
cepted negative c values as such. 

4.1 .2. Hubble diagram in UB only 

We applied the same procedure as in the previous section to fit 
the U and B-band light curves of the sub-sample of table[2]for 
which U -band measurements are available and redshifts larger 
than 0.015 (9 supernovae). We obtain M™ = -19.37 ± 0.05, 
a = 0.8 ± 0.4, B = 3.6 ± 0.6, and the standard deviation of 
residuals is 0.16±0.05. As expected, these results are consistent 
with the fit using B and V, as shown by the confidence contours 
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Fig. 5 68% joint confidence regions for (a, ft) fitted using either 
UB or BV light curves of the test sample. The crosses show the 
best fitted values with 1 <x uncertainties. 



for a and B fitted using either UB or BV light curves shown 
figure |5] Note the covariance between the estimated values of 
a and B, particularly in the U + B band case, which simply 
reflects the correlation between the parameters c and s in the 
test sample. 

One of the SN present in this sample (1996bo) show a large 
statistical uncertainty due to its limited [/-band photometry. 
Removing this point from the fit has the effect of bringing down 
the U, B contour to the point that it contains almost all of the 
B, V contour, bringing the two determinations of a and B closer 
to each other. 

Figure [6] presents the residuals to the Hubble diagram as 
a function of redshift, stretch and color using the values of 
M' B , a,B fitted with B- and V-band light curves in the previous 
section. Note the large uncertainty affecting 1996bo, which ap- 
pears as a ~ 2 sigma outlier in the 3 plots. Fitting the Hubble di- 
agram with the values of a and B obtained with B- and V-band 
light curves, the standard deviation of residuals is 0.20 + 0.05. 

Comparisons of the fitted values of s and c fitted using ei- 
ther UB or BV data are shown figure The error bars only 
reflect the propagation of photometric errors and do not ac- 
count for any intrinsic dispersion. There is no significant bias 
between the two estimates. 

4.2. Comparison with other luminosity distance 
estimators 

We compare our method with various estimators of SN la dis- 
tances by examining the dispersions about the Hubble line of a 
sample of nearby supernovae. The comparison is made on dis- 
tances measured in B and V bands only for which we obtain 
a dispersion of 0.16 ± 0.03. Our U- plus B-band distances es- 
timates could not be compared with other estimate due to the 
lack of published distances measured in these 2 bands. 

As shown in many papers (e.g. Wang et alJ 120031 and 
IWang et alJEooH) . testing luminosity distance indicators on 
"low-extinction sample" greatly improves their performance. 
As said in the introduction, cutting on color estimate measured 
with a redshift dependent accuracy is a source of systematic er- 
rors which we need to avoid for cosmological applications. We 
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Fig. 6 From top to bottom: Residuals to the Hubble diagram as 
a function of redshift, stretch s and color c indexes for super- 
novae of the test sample fitted in U and B bands. The values 
of M' B ,a,B used here are those fitted using B and V bands as 
described in the text. SNe with redshifts smaller than 0.015 are 
labeled with opened symbols. Note the large uncertainty affect- 
ing one SN ( 1 996bo) due to its limited [/-band photometry, and 
which also appears as a ~ 2 sigma outlier in the 3 plots. 



hence compare our resolutions to "full-sample" resolutions, for 
distance indicators involving rest-frame U, B and V bands. 

The MLCS method was originally presented in lRiess et alJ 
dl996ah . quoting a dist ance resolut ion of 0.12. Its latest devel- 
opment is presented in Jha ( 2002), with a distance resolution 
of 0.18, and used in Rie ss et al.l 112004). From the latter, we col- 
lected the distance measurements to the 20 objects in common 
with our test sample, and compute a Hubble diagram dispersion 
of 0.24 (0.22 with the low-redshift "golden sample") ) to be 
compared with our value of 0.16 when measured on the same 
20 events. 
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Fig. 7 Differences of 5 and c values fitted with either UB or 
light curves of supernovae from the training and the test 
sample (respectively opened and filled symbols). The error bars 
only reflect the propagation of photometric errors and do not 
account for any intrinsic dispersion. 



The CMAGIC method of IWang et alJ J2003h finds a 
weighted dispersion of 0.08 for a sub-sample of SNe with 
B max - V max < 0.05. With a weaker cut on color, B max - V max < 
0.5, the dispersion rises to about 0.15, which is co nsistent with 
our re sult. Similarly, the ACn method presented in lWang et alJ 
(2005), also present an exquisite distance resolution of 0.07 
in the V band based on a low-extinction subsample. However, 
when considering the full sample, the distance resolution de- 
grades to 0.18. A more detailed comparison with CMAGIC 
and AC12 would require comparing the distances of the same 
sample of SNe (due to the limited statistics) but their distances 
are not published. Note also that in these papers, the test and 
training samples are not separated. 

In summary, the method we propose gives a dispersion on 
distances measured using B- and V-band data only comparable 
or lower than obtained with other methods while also provid- 
ing, for the first time, comparable dispersion values for dis- 
tances measured using U- and B-band data only. 



high-redshift SNe la for which limited coverage is obtained in 
both wavelength and phase. 

The k-corrections, which allow the observer to transform 
the observed magnitudes into the standard rest-frame mag- 
nitudes, are built-in; the model includes the dependence on 
stretch and color of the spectrum template needed to estimate 
those corrections. In particular, the well-known correlation be- 
tween (U - B) max and stretch is reproduced. 

The (B - V) and (V - R) stretch-independent colors we ob- 
tain are extremely similar to the ones expected from reddening 
by dust. The (U - B) color departs from this law. We find a 
relation between (B - V) color and observed B luminosity in- 
compatible with Rb = 4.1, at more than 3 standard deviations, 
even when accounting for the stretch-color correlation. 

We have tested this fitting procedure on an independent 
sample of SNe la. Alternatively using B- and V-band data and 
U- and B-band data, we are able to retrieve consistent pa- 
rameters and hence build Hubble diagrams with both sets of 
data. The dispersions about the Hubble line were found to be 
0.16 + 0.03 and 0.16 + 0.05 in the B plus V and U plus B bands 
only, respectively. This method is particularly well adapted to 
reliably measure SN la distances in the full redshift range and 
in particular beyond redshift z ~ 0.8 for which rest-frame V- 
band measurements are often not available. 

Acknowledgements. It is always a pleasure to acknowledge the stim- 
ulating discussions among the FROGS (FRench Observing Group of 
Supernovae), especially with G. Garavini, J. Rich and R. Taillet whom 
we would also like to thank for their critical reading of the manuscript. 
We thank P. Nugent for providing us with his template spectrum time 
series. 



5. Conclusion 

We have proposed a new method to fit broadband light curves 
of Type la supernovae. It allows us to determine simultaneously 
the SN la rest-frame B magnitude at maximum, stretch and 
color excess (or deficit) using any measured multi-color light 
curve within the wavelength range of rest-frame UBV bands. 
This technique is particularly well-suited to the treatment of 
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Table 1 . The training sample of Type la supernovae light curves 



Name 


z a 


Bands 


m" B 


s 


c 


Phot. Ref. b 


1981B 


0.006 


UBV 


12.018 (0.003) 


0.911 (0.004) 


0.165 (0.002) 


(B83) 


1984A 


-0.001 


UBV 


12.389 (0.003) 


0.946 (0.005) 


0.215 (0.003) 


(K86) 


1986G 


0.002 


BV 


12.026 (0.006) 


0.736 (0.005) 


0.915 (0.006) 


(P87) 


1990N 


0.003 


UBVR 


12.701 (0.006) 


1.055 (0.005) 


0.090 (0.005) 


(L91,L98) 


1991T 


0.006 


UBVR 


11.574 (0.003) 


1.129 (0.005) 


0.183 (0.003) 


(F92,L98,A04,K04b) 


1992 A 


0.006 


UBVR 


12.546 (0.003) 


0.794 (0.003) 


0.088 (0.002) 


(S92,A04) 


1992al 


0.015 


BVR 


14.448 (0.012) 


0.922 (0.011) 


-0.035 (0.012) 


(H96) 


1994D 


0.001 


UBVR 


11.722 (0.002) 


0.780 (0.002) 


-0.068 (0.002) 


(R95,P96,M96,A04) 


1994S 


0.015 


BVR 


14.764 (0.017) 


1.006 (0.025) 


0.011 (0.018) 


(R99) 


1994ae 


0.004 


BVR 


13.108 (0.008) 


0.990 (0.006) 


0.098 (0.009) 


(R99A04) 


1995D 


0.007 


BVR 


13.450 (0.013) 


1.029 (0.014) 


0.072 (0.012) 


(R99,A04) 


1995al 


0.005 


BVR 


13.309 (0.016) 


1.038 (0.019) 


0.168 (0.017) 


(R99) 


1996X 


0.007 


BVR 


12.986 (0.010) 


0.868 (0.011) 


0.050 (0.010) 


(R99) 


1997E 


0.013 


UBVR 


15.101 (0.006) 


0.820 (0.010) 


0.078 (0.006) 


(J02) 


1997do 


0.010 


UBVR 


14.314 (0.014) 


0.920 (0.013) 


0.118 (0.009) 


(J02) 


1998bu 


0.003 


UBVR 


12.069 (0.002) 


0.989 (0.004) 


0.260 (0.002) 


(S99) 


1998dh 


0.009 


UBVR 


13.825 (0.011) 


0.861 (0.006) 


0.110 (0.008) 


(J02) 


1998es 


0.011 


UBVR 


13.814 (0.007) 


1.061 (0.009) 


0.104 (0.006) 


(J02) 


1999aa 


0.014 


UBVR 


14.728 (0.005) 


1.055 (0.005) 


-0.007 (0.005) 


(J02,A04,K00) 


1999ac 


0.009 


UBVR 


14.114 (0.005) 


0.925 (0.009) 


0.112 (0.006) 


(J02) 


1999cc 


0.031 


UBVR 


16.802 (0.009) 


0.834 (0.012) 


0.047 (0.010) 


(J02) 


1999cl 


0.008 


UBVR 


14.819 (0.008) 


0.901 (0.011) 


1.118 (0.007) 


(J02,K00) 


1999dq 


0.014 


UBVR 


14.398 (0.004) 


1.057 (0.006) 


0.118 (0.004) 


(J02) 


1999ee 


0.011 


UBVR 


14.839 (0.003) 


0.979 (0.003) 


0.297 (0.002) 


(S02) 


1999ek 


0.018 


UBVR 


15.616 (0.005) 


0.888 (0.007) 


0.167 (0.005) 
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"Heliocentric redshift. 

b Photometry Ref erenc es : B83: feuta & Turned ll983l). K86: iKimeridze & Tsvetkovl jl98d) P87: 
Phillips et alJll987l). L 91: lLeibundgut et alJll99ll> . F92: lFilippenko et alJll992l) . S 92:ISuntzefl]ll992l) . R95: 
Richmond et alJll995l) H96:lHamuv et alJll996ah. M96:lMeikle et allll99Cft P96:|Patat et alii 199(3). L98: 
Lira et alJ 119981) R99: [Riess et alJ ll999l). S99 : JSuntzeffet alJ Il999h. H01: iHowelll 1200 ll) L01: iLi et al. 
| 2001bl) . J02: IjhtJ 120021) S02: IStritzinser et alJ l2002h. K0 3: iKrisciunas et alj l2003h. V03: IValentini et 
| 2003 | ) . Vi03 :lvink6 et alJ 120031), Z03: IZapata et alJ 120031), A04: lAltavilla et alJ 12004]) BO^jTlBenetti et al. 
120041) . G04: lGaravini et alJl2004l) . K04b: IKrisciunas et alJl2004 . G05: lGaravini et alJl2005t) 
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Table 2. The test sample of Type la supernovae light curves 
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"Heliocentric redshift. 

b CMB-centric redshift, 
c Pho tometr y Reference s : H96: Hamuv^yd 



J2001bh. J02: Ijhl j2002j) . S02: IStrolger et al 
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2003), Z03: Zapata et al. 
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